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ABSTRACT 

We study the effect of massive neutrinos on the formation and evolution of 
the first filaments containing the first star-forming halos of mass M ~ lO^M© at 
z ~ 20. With the help of the extended Press-Schechter formalism, we evaluate an- 
alytically the rates of merging of the first star-forming halos into zero- dimensional 
larger halos and one-dimensional first filaments. It is shown that as the neutrino 
mass fraction fi, increases, the halo-to-filament merging rate increases while the 
halo-to-halo merging rate decreases sharply. For /j^ < 0.04, the halo-to-filament 
merging rate is negligibly low at all filament mass scales, while for fy > 0.07 
the halo-to- filament merging rate exceeds 0.1 at the characteristic filament mass 
scale of ~ IO^Mq. The distribution of the redshifts at which the first filaments 
ultimately collapse along their longest axes is derived and found to have a sharp 
maximum at 2; ~ 8. We also investigate the formation and evolution of the 
second generation filaments which contain the first galaxies of mass IO^Mq at 
2; = 8 as the parent of the first generation filaments. A similar trend is found: For 

> 0.07 the rate of clustering of the first galaxies into the second-generation fila- 
ments exceeds 0.3 at the characteristic mass scale of ~ lO^^M©. The longest-axis 
collapse of these second-generation filaments are found to occur at 2 ~ 3. The 
implications of our results on the formation of massive high- 2; galaxies and the 
early metal enrichment of the intergalactic media by supernova-driven outflows, 
and possibility of constraining the neutrino mass from the mass distribution of 
the high- 2; central blackholes are discussed. 



Subject headings: cosmology:theory — large scale structure of Universe 
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INTRODUCTION 



In 1970, Zel'dovich suggested that the large-scale structure of the Universe such as sheets 
and filaments form through the anisotropic collaps e of large-scale de nsity fluctuations along 
the principal axes of the gravitational tidal fields ( IZel'Dovichl Il970[ ) . Although his original 
model assumed a hot dark matter cosmology, this picture can be well acco mmodated even b y 
the currently popular A-dominated cold dark matter (ACDM) cosmology ( iBond et al.lll996l ). 
The dark matter halos that condense out first on the smallest scales merge anisotropically 
along the principal axes of the tidal fields into larger and larger structures. The inherent 
anisotropic nature of the gravitational merging process leads to the formation of filament- 
like and sheet-like structures, as envisioned by Zel'dovich's model. The nonlinear evolution 
of the tidal shear fields sharpens the anisotropic interconnection of the large-scale structure 
(which is often called the cosmic web), developing a hierarchical system w here small filaments 



are threaded with the larger parent filaments (e.g, iSpringel et al.ll2005l ). To understand the 
formation and evolution of the large-scale structure, an important question to answer is how 
frequently the halos merge into filaments and under what circumstances the halo-to-filament 
merging event occurs efficiently. 



Lacey fc Cold (119931 119941 ) have for the first time evaluated analyt ically the merging 



rates of bound halos by extending the original Press-Schechter formalism (jPress fc Schechtei 
19741 ). Their analytic work was later ref ined and complimente d by nume rical calculation 



of mer ger trees from N-body simulations (ISomerville et al.ll2000[ ). Recently, iFakhouri &: Ma 



( 20081) constructed merger trees from the large Millennium Run halo catalogs ( ISpringel et al. 



20051 ) and found that the mean mergi ng rate per halo can be well modeled by a nearly 
universal fitting formula. In addition, iFakhouri &: Mai (|2008[ ) showed that the prediction 
of the extended Press-Schechter formalism for the halo merging rates disagrees with the 
numerical results from the Millennium Run halo catalog by up to a factor of a few. 

The halo-to halo merging rates were originally derived under the assumption that the 
merging events occur in a hierarchical way. If the matter content of the Universe were 
composed purely of CDM particles, then the structure formation would proceed in a strictly 
hierarchical way. However, if the non-CDM particles (either warm or hot) coexist with the 
CDM particles, the pure hierarchy in the build-up of th e large-scale structu re would break 
down at some level, depending on their mass fraction ( lYoshida et al.l l2003l ) . True as it is 
that the current observations put a tight constraint on the amount of the possible non-CDM 
particles, they have not been completely excluded. As a matter of fact, the existence of 
non-CDM particles have been suggested as one possible solution to alleviate the apparent 
mismatche s between observations and the predictions of ACDM cosmology on the subgala ctic 
scales (see iBoyanovsky et al.ll2008l : iKuzio de Naray et al.ll2010t Ide Vega fc Sancheal2010l . for 
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a recent discussion). 



Among the various candidates for the non-CDM particles, the massive neutrinos have so 
far attained the most serious attentions. Ever since it was discovered that at least one flavor 
of the neutrinos must be massive and that the neutrino masses would be best constrained 



from the astronomical data (see lLesgourgues fc Pastoiil2006l . for a review), a flurry of research 
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In previ- 



ous studies of the structure formation in a AMDM (A+CDM+massive neutrinos) universe, 
what has been largely focused is the suppression of the small-scale powers in the density 
power spectrum due to massive neutrino's free-streaming and its effect on the structure 
formation. 

Since the massive neutrinos would decrease the merging rates of the small-scale halos due 
to their free streaming, one would naively think that the formation of large scale structures 
would be deferred in a AMDM universe. Meanwhile, the recent discovery of high- z quasars 



2009 



at ^ > 6 fjFan et al.l l200ll: lFanll2003l : iFan et al.l 12004 I2OO6I : I Jiang et al.ll2009l : IWillott et al. 



201of r has suggested that the first galaxies of mass > lO^M© must have formed at 



redshifts z > 8. At first glance the existence of such massive first galaxies at such early 
epochs might look difficult to reconcile with a AMDM model. Even in a ACDM cosmology, 
it is not so easy to explain the existence of the high-z quasars since it requires very high 
merging/accretion rate in the early Universe. 

Yet, what has not been noted in previous works is that the presence of massive neutrinos 
may actually open a new channel for the rapid structure formation. Since the massive 
neutrinos would behave as effective warm dark matter (WDM) particles, they would fasten 
the collapse of matter along the major axes of the gravitational tidal fields on large scales. In 
other words, while the free-streaming of massive neutrinos slows down the formation of the 
zero-dimensional bound halos on small scales, it speeds up the formation of one-dimensional 
filaments on larger scales. The filaments are the most optimal environments for the rapid 
accretion of gas and matter. Furthermore, when the filaments collapse along their longest 
axes, they would provoke the violent merging of the constituent halos. 

Here, assuming a AMDM universe, we study how the first filaments form through the 
clustering of the first star-forming halos and how they evolve after the longest-axis col- 
lapse. The most idealistic approach for the study of the first filaments would be to use 
high-resolution N-body simulations for a AMDM universe. But, it is still a daunting task 
to implement the dynamics of the massive neutrinos into the current N-body simulations. 
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Instead, we choose an analytic approach base d on the extended Press- Schechter formahsm 
( jPress fc Schechterl Il974j : iLacey fc Cold Il993l ) for this study. Throughout this paper, we 
assume fiat AMDM models where the neutrinos have mass and three species = 3. The 
key cosmological parameters are set at Q^n = 0.266, = 0.734, = 0.0449, ag = 0.801 
(for the ACDM part) , h = 0.710 and Us = 0.963, to be consistent with WMAP7 resuhs 
jKomatsu et al.ll201l[ l. 

The outline of this Paper is as follows. In §2, we derive analytically the rates of merging 
of the first star forming halos into the first filaments and determine their characteristic mass 
scales. In §3.1, we evaluate the epochs of the longest-axis collapse of the first filaments and 
the rates of merging of the first galaxies into the second-generation filaments, and determine 
their characteristic mass scales. In §4, we explore the possibility of constraining the neutrino 
mass from the observable mass distribution of the high-z supermassive blackholes. In §5, 
the results are summarized and a final conclusion is drawn. 



FORMATION OF THE FIRST FILAMENTS 



Recent observations indicate that the first galaxies must have forme d in the halos of 



mass 10^ - lO^M© at redshifts ^ = 9 - 10 (e.g., IWillott et al.ll2009l |2010|). What has been 
inferred from the current high-resolution hydrodynamic simulation is that these first galaxies 
may have f ormed through rapi d merging of the first-star forming minihalos of mass lO^M© 
at z > 20 (IBromm et al.ll2009l . and references therein). In this section, we investigate how 
the presence of massive neutrino would change the evolution channel of the first-star forming 
halos. 



2.1. Variation of the Halo-to-Halo Merging Rates with the Neutrino Mass 



In the extended Press-Schechter theory developed by iLacey fc Cold (Il993l . Il994j ). the 
key quantity is the conditional probability, f{M/, z/\M, z), that a halo of mass M existing at 
a given redshift z will merge to form a more massive halo of mass M/(> M) at some lower 
redshift z/{< z). Adopting the spherical collapse condition that a halo of mass M forms at z 
when its linear density contrast 6 on the mass scale M reaches a redshift-depe n dent threshold 
5c/D[z) where 5c ~ 1.68 and D{z) is the linear growth factor. iLacey &: Cold (119931 ) equated 
this key conditional probability to the differential fraction of the volume in the linear density 
field occupied by the regions satisfying 5i > 6c/D{z/) on the mass scale M/ at z/ provided 
that it also satisfies the condition of 5 = 6c/D{z) on the mass scale M at earlier epoch z. 
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This key conditional probability f{M/, z/\M, z) dMi can now be written as 



f{S!,u!\S,u)dSt 



2^p(5a/^ > uj'\5m = w) dSf 



dSi 



v[.^Mi\^M = uj)d6M,dSf 



S 



St{S - St) 



exp 



{ujiS - uSif 
'2SSi{S-Si) 



dSf, (1) 



where /(M/, zf\M, z) dMi = f{Sf, ujf\S, oj) dSi with S = a'^{M),u = 5,(2), Si = a^iMf), and 
u/ = 6c{z/). Here a{M) and a{M/) are the rms fluctuation of the linear density field smoothed 
on the mass scale M and M/, respectively. Normalization of a is done for the AMDM linear 
power spectrum to be identical with the ACDM linear power spectrum at large scale. Here, 
of course, with the normalized ACDM linear power spectrum a at 8h~^Mpc has the value of 
0-8 given by WMAP7. 



Using Equation (fB. lLacey &: Cold (Il993l ) calculated the (instantaneous) rate of merging 
of a dark halo with mass M at a given redshift z into a larger halo with mass M/ by taking 
the derivative of f{Sf,ijjf\S,ijj) with respect to zf: 



d^M Ml\z) 
dlnAMdz 



dSf 
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dz/ duji 
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3/2 



Sl{S - Sf) 



exp 



1 



(2) 
,(3) 



where AM = Ml—M and -Pl(/c) is the linear d ensity power spect r um fo r a AMDM cosmology. 
We use the analytic approximation given by IHu fc EisensteinI (119981 ) for the evaluation of 

PL{k). 

With the help of the above analytic prescriptions, we determine the effect of massive 
neutrinos on the merging rate of the first star-forming mini-halos of mass M = 10^ Mq 
at redshift z = 20 for the ra nge of the neutrino mas s fraction, 0.0 < fy < 0.13. We use 
the analytic formula given by iHu fc EisensteinI ( 1l998l ) for the AMDM linear density power 
spectrum = Q^/Qm where Qi, is the neutrino mass density parameter. The AMDM linear 
growth factor D{z) is appro ximated as D^np^^^f^r) wher e -Dacdm(^) is the linear growth 
factor for a ACDM universe (ILesgourgues fc Pastoiil2006l ). The range of the neutrino mass 
fraction is chosen to be 0.0 < /„ < 0.13 wh ere the upper limit of 0.13 comes from the recent 
WMAP7 constraint ( iKomatsu et al.l 120111 ). It is also worth mentioning here that for this 
range of fi, the mini-halos of mass IO^Mq at z = 20 correspond to the (3 — 4)a peaks of the 
initial density field. 
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The rates of the merging of the first star-forming mini-halos are plotted as thin hnes 
in Figure [1] for the five different cases of f,^, which reveals that the halo-to-halo merging 
rates decrease rapidly as fi, increases. This result implies that in a AMDM universe the free 
streaming of the massive neutrinos have an effect of deferring the formation of larger halos 
through the major mergers of the first star-forming mini-halos. 



2.2. Variation of the Halo-to-Filament Merging Rate with the Neutrino Mass 



The result of §2 indicate that the direct merging of the first mini-halos into larger halos 
are slowed down by the free streaming of massive neutrinos. But, what about the merging of 
the first star-forming mini-halos into the first mini-filaments? To derive the rate of merging 
of the first st ar-forming rn i ni-ha los into the mini-filaments, we adopt the classification scheme 
proposed by lForero et al.l (120091 ). according to whom a filament of mass M/ forms at z/ when 
the two of the three eigenvalues, A/i, A/2, A/3, of the linear deformation tensor (defined as 
the second derivative of the linear peculiar potential) smoothed on the scale of M/ exceed the 
threshold value, Xc/D{zf) with Ac ~ 0.3. Using the N-body results, they demonstrated that 
the structures identified by this dynamical classification scheme exhibit indeed filamentary 
sh apes. Their filament-find ing scheme is also consistent with the statistical results obtained 
by lLee &: ShandarinI (119981 ). 



Adopting this condition for the formation of filaments, we modify the Lacey-Cole for- 
malism to determine the halo-to-filament merging rates. In this modified formalism, the key 
quantity is the conditional probability, /(M/, z/\M, z), that a halo of mass M observed at z 
will merge into a filament of mass M/ formed at z/. We equate this conditional probability 
to the differential volume fraction occupied by the regions satisfying the filament formation 
condition, A/i > A/2 > Xc/D{zf), A/3 < Xc/D{zf) on the mass scale Mf at redshift zf provided 
that the linear density contrast of the same regions, 6, on the smaller mass scale M at earlier 
epoch z satisfies the spherical collapse condition 6 = 6c/D{z). Here 6 equals the sum of the 
three eigenvalues of the deformation tensor smoothed on the mass scale M. 



Now, this key conditional probability can be expressed as 

f(Sf, r]f\S, oj) dSi = 2-^p(A/i, A/2 > 77/, A/3 < r]f\5 = co) dSf, 
dSf 



(4) 



with rjf = Xc/D{zl). To evaluate f{Sf,rif\S,u) dSf in Equation (jl]), it requires us the joint 
conditiona l probability distribution ofpfAi, A2, A 3 1^), w hich has been already derived by 
Led ( 120061 ) from the original work of iDoroshkevichl (jl970[ ) as 

p(A/i,A/2,A/3,5) 



p(A/i,A/2,A/3|5) 



p{6) 
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A/2 + A/3 and II2 = 



(5) 



A/1A/2 + 



It is now straightforward to calculate the (instantaneous) rate of merging of the first 
star forming halos of mass M = 10^ Mq at 2; = 20 into the early mini-filaments of mass M/ 
as 



d^pjM Mt\z) 
rflnAMrfz 



AM ^ . d 



dzi 



(6) 



Figure [Uplots the halo-to-filament merging rates (thick lines). In contrast to the halo-to-halo 
merging rates, the halo-to-filament merging rates increase with f^. As fy changes from to 
0.13, the maximum halo-to-filament merging rates increase by a factor of 3. We call the scale 
on which the halo-to-filament merging rate reaches the maximum value the characteristic 
filament mass scale and denote it as Mp- For > 0.07, the halo-to- filament merging rate 
begins to exceed 0.1 at the characteristic mass scale around Mp = 10^ Mq. Whereas for 
fi, < 0.04, the filament-to-halo merging rate is below 0.1 at all filament mass scale. Note 
that for fi, = 0.07 the mini-halos of mass 10^ Mq at 2; = 20 tend to merge more rapidly into 
the mini- filaments of mass 10^ Mq rather than into the larger halos of mass 10'' Mq. Figure 
|2] plots the characteristic mass Mp of the first filaments versus the neutrino mass fraction 
/j,. As one can see, the value of M^, decreases almost linearly with f^. 



3. EVOLUTION OF THE FIRST FILAMENTS 

3.1. Longest-Axis Collapse of the First Filaments 

In the previous section, it is shown that the first star-forming mini-halos would merger 
rapidly into the first filament of characteristic mass ~ W^Mq in the presence of massive 
neutrinos. As the Universe evolves, these first filaments would ultimately collapse along 
their longest-axes that are in the direction of the minor principal axes of the local tidal 
field. To determine the distribution of the redshift Zf when the first filaments collapse along 
their longest axes, we calculate the conditional probability density that a region satisfying 
the filament-formation condition (Ai > r], A2 = rj, < rj) at z = 20 on the mass scale Mp 



- 8 - 



will meet the halo-formation condition {5 = u) at the same mass scale but at some lower 
redshift, zj <20: 



p{Sf = Uf\\i >r],\2 = T], As < T]) 



I-oo ^-^1 P(^' = -^1' -^2 = V, ^3 

Hoc XT ^2 = V, h) 



(7) 



where S/ is the linear density contrast on the mass scale Mp at Zf, Ai, A2, A3 are the three 
eigenvalues of the deformation tensor smoothed on the mass scale of Mp aX z = 20, ui = 
Sc/D{zf) a nd r? = XrJD(z ) . Usi n g the loint p robability density distribution of p{Sf, Ai, A2, A3) 
derived by iDoroshkevichI (119701 ) ; iLed (120061 ) , it is also straightforward to evaluate equation 



Figure [3] plots this distribution p{zf) of the epochs of the longest axis collapse of the 
first filaments for the five different cases of f^. As can be seen, p{zf) has a maximum value 
at Zf = 8-9 for all five cases, which implies that the first filaments retain their structures 
for a long period from ^ = 20 to 2; = 8-9 without collapsing into zero- dimensional halos. As 
/jy increases, however, the distribution p{zf) becomes narrower and its peak position moves 
slightly to the lower value of 2;/. That is, the massive neutrinos play a role of delaying the 
longest-axis collapse of the first filaments. 



3.2. Formation and Evolution of the Second-Generation Filaments 

Now that the halo-to-filament merging rates are found to increase with the neutrino 
mass fraction, we would like to see how the presence of the massive neutrinos alters the 
merging rates of the first galaxies. It is naturally expected that some fraction of the first 
galaxies would also merge first in the filaments (called the second-generation filaments). 
A question is what the characteristic masses of the second generation filaments are and 
when they would undergo the longest-axis collapse. To answer this question, we recalculate 
everything through equations [UlTl but for the first galactic halos of mass 10^ Mq at 2; = 8. 

Figure H] plots the halo-to-halo and halo-to-filament merging rates of the first galactic 
halos of mass 10^ Mq at 2 = 8 as thin and thick lines, respectively. A similar trend to 
that shown in Figure [1] is found for the first galactic halos: As the neutrino mass fraction 
increases, the halo-to-halo merging rate decreases while the halo-to-filament merging rate 
increases. The comparison of Figs{T] and H] also reveals that the halo-to-filament merging rates 
increases with redshift and with halo mass. For > 0.04, the rate of merging of the first 



To prevent the mathematical divergence, the mass of a final halo is set at 0.99Mf. 
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galaxies into the second generation filaments exceeds 0.1 at all mass scales. The characteristic 
mass of the second generation filaments is found to lie in range of 10^^ < M/Mq < 10^^. 

Figure E] plots the distribution of the redshifts at which the second generation filaments 
of mass IO^^Mq formed at z = 8 experience the longest axis collapse. A similar trend to 
that shown in Figure [3] is also found: p{zf) has a maximum value at Zf = 3-4 for all five 
cases, which implies that the second-generation filaments last from z = 8 to z = 3-4 without 
collapsing into zero-dimensional halos. The increase of fi, results in narrowing down the 
shape of p{zf) and shifting it slightly to the low-z section. Note that these second generation 
filaments correspond to filaments hosting galactic halos and are ex pected to be interconn ected 



to the larger-scale parent filaments of mass > 10^^ — IO^^Mq (e.g.. ISpringel et al.ll2005l ). The 
massive galaxies formed through the longest-axis collapse of the second generation filaments 
would grow through the anisotropic accretion of matter and gas along the larger-scale parent 
filaments in the cosmic web. 



CONSTRAINING THE NEUTRINO MASS WITH THE FILAMENT 

ABUNDANCE 



Once the first filaments form through the merging of the first star-forming halos at 
z < 20, the first stars comprising the first filaments would grow more rapidly as t he accretion 
of ma tter and gas occur more efficiently along the bridges of the first filaments ( iPark fc Lee 
20091 ). When the first filaments collapse along their longest axes at 2; ~ 8, the strong collisions 
between stars and gas particles in the filaments would result in seeding the supermassive 
blackhol es in the massive fir st galaxies that are believed to power the ultra-luminous high- 2; 
quasars ( IWillott et al.ll2010l ). In this scenario, the high-end slope of the mass distribution 
of the supermassive blackholes inferred from the observable luminosity function of the ultra- 
luminous high-2 quasars would refiect the mass distribution of the first filaments. Since the 
mass distribution of the first filaments depend sensitively on the mass of massive neutrinos, 
it might be possible to constrain the mass of neutrinos by comparing the mass distributions 
of the first filaments with that of the high- 2; supermassive blackholes. 

The mass distribution function of the first filaments {dNp/dM) at the epoch of longest- 
axis c ollapse (zf) can be derived by employing the Press-Schechter approach (jPress fc Schechter 
193): 



dNpiM, Zf] 
dM 



M 



d 
dM 



F{M, Zf) 



(8) 



where p is the mean mass density. Here, F{M, Zf) represents the volume fraction occupied 
by the first filaments of mass larger than M at Zf. Following the Press-Schechter approach. 
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this volume fraction can be expressed as 

F{M,Zf) = Ap[X^ > A2 > Ae,.,,A3 < Ae,.>M], (9) 

where A is the normahzation factor, A^ is the value of the A2 at the epoch of the longest- 
axis collapse. Since the first filaments formed at the threshold of Ac = 0.3, the value of A2 
increases gradually till they collapse eventually along their longest axes. It has been found 
that at the epoch of the longest-axis collapse, A2 reaches the critical value of A^ = 1.0. 



Using the joint distribution, of p(Ai, A2, A3), derived by iDoroshkevichl (Il970l ). it is 
straightforward to evaluate Equation ([8]). Figure [6] plots the mass distribution of the first 
filaments in logarithmic scales at the epoch of the longes-axis collapse for five different cases 
of the neutrino mass fraction, f^. As we have already shown in §2, that the longest-axis 
collapse of the first filaments occur around z = 8, we set at 8 for the evaluation of 
Equation ([8]). As it can be seen, in the high-mass section (M > 10^ /i^^M©), the abundance 
of the first filaments decreases sharply with M for all cases. The rate of the decrease of the 
mass function of the first filaments, however, depends sensitively on the value of /j,. The 
larger the value of f^, is, the more rapidly the mass function drops in the high-mass section. 

To examine quantitatively how the rate of the decrease of the mass function of the first 
filaments in the high- mass section changes with fi^, we approximate it as a power law and 
determine the value of the high-end slope, a, in the mass range of 10^ < M/ [h'^Mo) < 10^^. 
Figure [7] plots the high-end slope, a, of the mass distribution of the first filaments versus the 
neutrino mass fraction, fy as solid line. As it can be seen, the absolute value of the high-end 
slope of the mass function of the first filaments increases sharply as fp increases. This result 
suggests that the value of fy could be constrained by using the high-end slope of the mass 
function of the first filaments. In practice, of course, it is not possible to measure directly 
the high-end slope of the mass function of the first filaments. As mentioned above, however, 
we can infer it from the observable high-z blackhole mass function, under the assumption 
that the longest-axis collapse of the first filaments with mass larger than the characteristic 
mass (10^ H'^Mq) would result in the first galaxies hosting the supermassive blackholes that 
power the ultra-luminous high- 2; quasars. 

Now, we would like to compare the high-end slope of this theoretically derived mass 
distribution of the first filaments with that of the observationally obtained ma ss function of 



the m assive galaxies hosting the high- 2; supermassive blackholes. Very recently, IWillott et al. 



(120101 ) have derived the blackhole mass function at ^ ~ 6 from the observed luminosity 
function of the ultra- luminous high-z quasars (see Figure 8 in Willott et al. 2010). Converting 
the mass of the high-z supermassive blackholes into the t otal mass of their host galaxies 



through the following relation given by lBandara et al.l (120091 ) . log(Mbh/M0) = (8.18±0.11) + 



1^1.55 ± 0.31) X (log(Mtot/M0) — 13.0), we find the best-fit high-end slope is in the range of 
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— 6 < a < —4.3, which are plotted as two dashed hnes in Figure [71 Comparing this range 
of a against the theoretical curve plotted in Figure [71 it is found that the neutrino mass 
fraction is constrained as 0.09 < ]v < 0.13 which corresponds to the neutrino mass range of 
0.375 < m^l{eV) < 0.542. 



5. DISCUSSION AND CONCLUSIONS 

Starting with the first star-forming halos of mass IO^Mq at z = 20, we have shown ana- 
lytically that the presence of massive neutrinos opens a new channel for the rapid formation 
of massive structures by speeding up the halo-to-filament merging process. Once the first fil- 
aments of mass < 10^ Mq form through the clustering of some of the first star-forming halos 
at z > 9, matter and gas would accrete more efficiently along the first filaments, resulting in 
the fast growth of the constituent halos in the first filaments. As the Universe evolves and 
the tidal forces increases, the first filaments would undergo the ultimate collapse into the 
zero- dimensional halos (corresponding to the first massive galaxies) in the directions parallel 
to the minor principal axes of the local tidal field. 

What has been found here is that the longest-axis collapses of the first filaments end up 
forming the first massive galaxies of mass < 1O^M0 at z > 8. The longest axis collapses of 
the first filaments would accompany violent mergers among the constituent halos, which in 
turn would trigger the vigorous star-formation and feed the central blackholes in the resulting 
first galaxies. Our result suggests that the ultraluminous high-2; quasars detected at z > 6 
might correspond to the first massive galaxies formed through the longest-axis collapses of 
the first filaments. 

It has been also found that the second-generation filaments of mass lO^M© form through 
the clustering of some of the first galaxies of mass lO^M© at z ~ 8. The first galaxies 
in the second ge neration filaments would evolve fast through the filamentary accretion of 



matter and gas ([Park fc Leell2009l ). When the second generation filaments collapse along 
their longest axes, they would also induce violent mergers among the constituent galaxies, 
which would lead to the formation of massive galaxies of mass larger than IO^^Mq with high 
star- formation rates and large central blackholes at redshifts as high as z = 3. 

We test the possibility of constraining the neutrino mass, m^, from the mass function of 
the high- 2; supermassive blackholes that power the ultraluminous high-^ quasars, under the 
assumption that the longest axis collapses of the first filaments would seed the supermassive 
blackholes. Comparing the high-end slope of the theoretically derived mass function of the 
first filaments at the epoch of their longest-axis collapse with that of the observed mass func- 
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tion of the massive first galaxies inferred from the luminosity function of the ultra-luminous 
high- 2; quasars, we find that the neutrino mass is constrained to be 0.375 < {niy/eV) < 0.542. 

Yet, it is still only a feasibility study. A more comprehensive work should be done to 
put a robust constraint on my from observational data. First of all, more refined theoretical 
models for the halo-to-filament merging rates and the abundance of the first filaments should 
be constructed than the crude approximations based on the extended Press-Schechter formal- 
ism, since the formalism has bee n shown to fail in acc urately predicting the merging rates 
and abundance of bound halos (IFakhouri fc Mall2008l ). Second, the observational results 
from the luminosity function of the ultra-luminous high-z quasars still suffer from small- 
nu mber statistics a nd th us more quasar data are required. Third, we use the relation given 
by iBandara et al.l (120091) t o cori vert the blackhole mass into the halo mass. However, the 
relation of iBandara et al.l (|2009[ ) has been obtained from the blackholes detected at z < 4, 
and thus it may not be applicable to the higher-2; blackholes. A more accurate relation 
between the high- 2; blackhole mass and host halo mass should be found. 

Our results also hints that the massive neutrinos ma y have something to do with the 
early metal enrichment of the intergalactic medium (IGM). lMadau et al.l ( 1200 ll ) showed that 
the early outflows driven by supernovae (SN) ejecta from the galaxies of mass > lO^M© 
at z ~ 9 would enrich the IGM with the products of stellar nucleosynthesis provided that 
those galaxies possess relatively high star-formation efficiencies. It is interesting to note that 
the halo mass and epoch required for the SN-driven outflows are more or less coincident 
with our estimates for the galactic halos formed through the longest-axis collapse of the first 
filaments in a AMDM model. Assuming that the longest-axis collapse of the first filaments 
trigger star formation rapidly in the first massive galaxies in a AMDM model, the star- 
formation efficiency of the first mass galaxies may be as high as required for the early metal 
enrichment by SN-driven outflows. Furthermore, given our speculation that the longest- 
axis collapse of the first filaments would feed the supermassive blackholes, the anisotropic 
outflows from the active galactic nuclei (AGN) powered by the supermassive blackholes of 
the fi rst massive galaxies might also contrib ute to the early metal enrichment of the IGM 



[e.g., iGermain et al.ll2009l : iBarai et al.l 1201 ll ). It will be intriguing to explore how the metal 
enrichment of the IGM depends on the neutrino mass. Our future work is in this direction. 
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Fig. 1. — Rates of the merging of the first star-forming halos with M — 10^ Mq into the 
first filaments (thick hnes) and larger halos (thin hnes) at 2; = 20 for the five different values 
of the neutrino mass fraction: fi, — 0.13, 0.10, 0.07, 0.04 and 0.0 (sohd, dotted, dashed, 
dot-dashed, and dot-dot-dashed, respectively). 
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Fig. 2. — Histogram for the characteristic masses of the first filaments versus the neutrino 
mass fraction /j,. 



- 18 - 



0.6 




z 



Fig. 3. — Probability density distribution of the redshifts when the longest-axis collapse 
of the first filaments occur for the first different cases of the neutrino mass fraction: 
fi, — 0.13, 0.10, 0.07, 0.04 and 0.0 (solid, dotted, dashed, dot-dashed, and dot-dot-dashed, 
respectively) . 
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Fig. 4. — Same as Figured] but for the first galactic halos of mass M = lO^M© at 2; = 8. 
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Fig. 5. — Same as Figure [3] but for the second-generation filaments of mass M = IO^^Mq. 
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Fig. 6. — Mass distribution of the first filaments for five different cases of the neutrino mass 
fraction /j,. 
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Fig. 7. — High-end slope of the mass function of the first filaments versus f,y at Zf = 8 
(solid line). Ob servational co i istrain ts from the mass function of the high-z supermassive 
blackholes from IWillott et al.l (|2010[ ) are shown as dashed lines. 



